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ABSTRACT 

We study QCD correction to the pair production of colored scalar particles 
in electron-positron annihilation with an emphasis on gluon emission in the final 
state. We discuss the usefulness of working in a "quasi-two-body" frame and 
present the helicity amplitudes for the process. We compare the final state con- 
figuration with fermion pair production and find that the three-jet fraction for 
the scalars shows quantitative difference from that for fermions. 



1. Introduction 



The known elementary particles are either spin-1/2 fermions (quarks and 
leptons) or spin-1 gauge bosons. No elementary spin-0 particle has been found 
to date. Yet, scalar particles appear in most field theory models of particles. 
Higgs bosons are the key ingredient for electroweak symmetry breaking in the 
standard model and its extensions. Supersymmetry predicts a scalar partner 
for each fermion degree of freedom, thus the existence of three generations of 
squarks and sleptons. Grand unified Eq model contains colored scalar particles 
which may be interpreted either as leptoquarks or diquarks. 

Many of these scalars can be pair produced in e + e~ annihilation with a cross 
section comparable to or somewhat smaller than that for fermions. Experimental 
searches for these particles have been extensively performed [1]. 

If the scalar particle is colored (like the squark and leptoquark), the cross 
section is modified by QCD corrections. The calculation of the 0(a s ) correction 
is essentially the same as scalar QED one-loop calculation which is known for a 
long time. 1 The 0(a s ) correction for the scalar pair production is numerically 
quite important. In the high energy limit, the total correction factor for scalar 
pair production is four times larger than that for fermion pair production [3]: 



where cxo is the lowest-order cross section and Cr is the second-order SU(3) 
Casimir eigenvalue (4/3 for the fundamental representation like squarks or lep- 
toquarks). We use the notation ( to represent a generic colored scalar particle in 
the representation R. At lower energies, especially in the threshold region, the 
correction factor becomes even larger. 

The total 0(a s ) correction consists of two parts, one-loop virtual correction 
to the lowest order process e + e~ — > CC an d real gluon emission correction e + e~ — > 

1 A detailed description can be found in Ref. 2. A typographical error is corrected in Ref. 3. 
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CCg, both of which are infrared divergent. This divergence is cancelled when the 
two contributions are added. Essential for this cancellation is only the part of 
the three-body phase space where the gluon is soft (or collinear to one of the 
scalars when the scalar mass goes to zero). The rest of the three-body phase 
space corresponds to the final state with a hard gluon which may be observed as 
a separate jet ("three-jet" final state). 

The gluon emission process in scalar top quark pair production has been 
studied by Beenakker, Hopker, and Zerwas [4], who calculated the fully differ- 
ential cross section for e + e~ — > tig. In the present paper, we derive the cross 
section in a different Lorentz frame ("quasi- two-body" frame) in which the he- 
licity amplitudes have a simple interpretation. We then extend their analysis 
and compare the three-jet cross sections with those for scalar and fermion pair 
production processes. 

The QCD correction we calculate in this paper is just the effects arising 
from gluon gauge interactions. In the supersymmetric standard model there are 
additional interactions with the same strength involving gluinos. The gluino 
coupling contribution to the virtual 0(a s ) correction is discussed by Arhrib, 
Capdequi-Peyranere, and Djouadi [5]. 

2. Lowest-Order Cross Section 

We consider the process e + e~ — > CC which occurs via s-channel 7 or Z ex- 
change. Additional t or w-channel exchanges have to be included for selectron 
(or electron sneutrino) production or a leptoquark if it couples to electrons with 
a nonnegligible Yukawa coupling. 

In this section we list the lowest-order amplitude and cross section for com- 
pleteness. The Feynman graph is shown in Fig. 1. The amplitude is found to be 
(see Fig. 1 for the momentum assignments) 

Mo = -H^(p-p)^, (2.1) 
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with 

H" = -Q C v(Pe)l»<Pe) + T 3C - Q < Bm *° W -^v(p e mv e ~ a el5 )u(p e ) . 

cos z 6\y sin d w s - 

(2.2) 

Here s is the e + e~ cm. energy squared, and are the electric charge and 
the third component of weak isospin of £, and 

v e = ~ + sin 2 e w , a e = ~. (2.3) 

We will neglect the electron mass throughout the paper. 

For longitudinally polarized electrons with polarization P colliding with un- 
polarized positrons, the cross section can be written 

1 _|_ p 1-P 

da = — - — da+ H — da- , (2.4) 

where da± denotes the cross section for the e^e£ and e^e^, initial states. We 
find 

da± nd R a 2 /3 3 2 . 2 



dcos9 2s 

with 



-Hi sin^ , (2.5) 



;W e =Fa e )(T 3C -(2 c sin 2 (%) s 



cos z 6w sin dw s ~ m z 
Here, d R is the dimension of the color SU(3) representation of C (3 for the fun- 
damental representation), and 



13 = yj\ -4m 2 /s , (2.7) 

with m being the ( mass. The f3 3 factor reflects the P wave threshold. The total 
cross section is 

2nd R a 2 /3 3 2 
v± = ^ H ± » ( 2 - 8 J 

which we will denote by ao±. 
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3. 0(ct s ) correction 

3.1 Virtual one-loop correction 

Feynman diagrams for the process e + e~ — > CC at 0(«s) are depicted in 
Fig. 2. The diagrams (a)-(c) are the one-particle-irreducible vertex corrections. 
The mixed four-point vertex in (b) and (c) appears because the scalar kinetic 
term has two covariant derivatives. The diagrams (d) and (e) are £ self-energy 
corrections, and the diagrams (f)-(h) show the counterterm contribution. 

2 

We regularize the ultraviolet divergence by dimensional reduction with D = 
4 — 2e, and the infrared singularity by an infinitesimal gluon mass A. This pro- 
cedure does not cause problems with gauge invariance since the diagrams do not 
contain non-Abelian gauge vertices. 

We adopt the on-shell renormalization scheme to determine the counterterms. 
The mass and wave function renormalization constants for ( are found to be 



5m 2 Cjia t 



m 2 Air 
Z c - 1 



1 , A „ , rn 



■2 



3( - -7i5 + log47rJ -31og-^- + 7 

1 \ A 2 

2 ( - - 7# + log47rJ - 2 log 



(3.1) 
(3.2) 



An 

where fx is the renormalization scale and is the Euler constant. 
We write the £C7 vertex function in the form —ieQ^A^ with 

A„ = F(q 2 )( P -p)^. (3.3) 

(See Fig. 2 for the definition of the momenta.) Up to the order we work, the ((Z 
vertex can be writtten as 

-ie(T K - Q c sin 2 ^) 

COS Vw Sin Vw 

with the same A^. At the lowest order we have F(q 2 ) = 1. 



2 Naive dimensional regularization gives the same physical results, although the finite parts 
of the renormalization constants differ. 
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The counterterm for the vertex can be found using the Ward identity. The 
same result is obtained by requiring no 0{a s ) correction at q 2 = 0, i.e., F(0) = 1. 
Expanding the form factor as 



F(q 2 ) = 1 + ^7(4 2 ) , (3-5) 



we find for q 2 > Am 2 



+ 



1+/3 

(3 



l-/3\ . 2/3 . 1+/3 1 . 2 1+/3 , 7T 



21 



Lio — log log log h 

1 1+0 ) S 1+/3 S 1-/3 4 s 1-/3 3 



1+/3 

+ — l ° g —p- 2 

1+/3 2 / m 2 , A/3 2 \ 
+ -^f (log ^ + log ^-2]. (3.6) 

3 x 

Here Li2(x) is the dilogarithm (Spence) function Li2(x) = — J Q dt log(l — t)/t. 
At 0(a s ), the lowest order cross section is multiplied by the factor 

l + ^3fe/(a). (3.7) 

7T 

The infrared singularity is cancelled by real gluon emission, to which we now 
turn. 



3 For a convenient expansion for numerical evaluation of dilogarithm, see Ref. 6. 
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3.2 Real gluon emission 



The Feynman diagrams for the process e + e — > ((g at the lowest order are 
shown in Fig. 3. The amplitudes are found to be 



{p-p + k)^p a _ (p-p- k)^p a 
p-k p-k 



~*9 t 



e* a . (3.8) 



Here T a is the SU(3) generator in the representation R. 

We calculate the cross section for this process in two methods, with agreeing 
results. One is the conventional trace technique. The other method is a helicity 
amplitude technique described in Appendix A. 

The differential cross section may be parametrized by four variables. We find 
it convenient to use the following four: 

t : CC c - m - energy squared normalized to s, r = (p + p) 2 / s; 

: angle between the e~ and the gluon in the e + e~ cm. frame; 

9 : angle between the £ and the gluon in the cm. frame; 

(j) : azimuthal angle between the e~ and the £ with respect to the gluon, 
measured from the e~ to the ( direction (common to both frames). 

The use of these variables is motivated by the fact that the amplitude can 
be split into two subprocesses e + e~ — > V* (V* = 7, Z) and V* — > ((g, and 
the latter process can be most conveniently evaluated in the CC c - m - frame. The 
detail and the result for the helicity amplitudes may be found in Appendix A. 

The polarized cross section is 
da± d R a 2 Hl C R a s 

V(l — T) 



dr d cos 6 d cos 9d(p 8s n 

x |a(1 + cos 2 6) +B(1-3cos 2 0) + Csin cos cos <j> + D sin 2 6 cos 20 j , 

(3.9) 



7 



where 



2/3 2 tv 2 sm 2 

(l-r) 2 (l-<; 2 cos 2 0) 2 ' ( } 

2r?; 4 sin 2 flcos 2 fl 

4y/rv 2 sin cos 9 [/3 2 - t> 2 + (1 + t)v 2 sin 2 6>] 
° = (l-r) 2 (l-^ 2 cos 2 ^) 2 ' (3 ' 10C) 

2rt; 2 sin 2 ^(/3 2 -^ 2 cos 2 ^) 

u - (i-T) 2 (i-v 2 cos 2 e) 2 • ^ iUdJ 



Here 



v = ^l-4m 2 /rs (3.11) 

is the velocity of ( in the (( cm. frame. 

If one integrates the cross section over the whole 3-body phase space, one 
encounters infrared divergence coming from soft-gluon region. We separate this 
region from the rest by the condition that the gluon energy k° < u, where the 
cutoff u is infinitesimally small. 

In this region, we regularize the divergence by giving an infinitesimal mass 
A to the gluon, which satisfies A <C u. We recalculate the amplitude with finite 
mass. The amplitude factorizes into the lowest-order part and the soft gluon 
factor. Since a very soft gluon does not alter the kinematics of the (( state, we can 
also integrate over the soft gluon phase space ignoring momentum conservation. 
Integrate over this soft-gluon region, we find the cross section 



C R a s \(l+0 2 ^ 1+0 \ 4u 2 

C-soft = CTQ < —— - log - 1 log ■ 



7T 1 V 2/3 to 1-13 / A 2 
l+(3 2 



+ 







fl-0\ 20 1+0 1 2 l+0 n 2 

U2 {lTi3)- l0g lTi3 l0g iq3-4 l0g 1=0 -T 
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1 , 1+P 



(3.12) 



where ao is the lowest-order cross section. 

For the rest of the phase space, we can set the gluon mass to zero. We first 
integrate over O, 4> to find 



do- 



dr d cos 6 



= O"0± 



C R a s v(l - t) 



n 



P Z 



1 + 



2(3 2 tv 2 sin 2 # 



(l-r) 2 (l-u 5 



cos 2 #)^ 



(3.13) 



which we will use to study the event configuration in the next section. The 
integrated cross section is 



C R a s \(l + p 2 . 1+P 



1 log 



m 
Aj 2 



+ 



1 + P 2 

P 



2Li2 (l^) +2Li2 



1-P \ 
1+P 



l+p 



l+p °i-p 



1 2 1+P n 2 

+ 2 log T^-¥ 



-31o gT -^ -41og/3 



4/3' 



3 (3 + P 2 ) (1 - P 2 ) log ^ + ^ (3 + 7/3 2 ) 



(3.14) 



The oj dependence cancels out when these two cross sections are summed. 
The infrared divergence is cancelled by adding the virtual one-loop correction 
discussed in the previous subsection. 
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3.3 Total correction 



The 0(a s ) corrected cross section for e + e — > (((g) can thus be written as 



a = ao 



1 + ^A 



TV 



(3.15) 



where 



1 



1 



A = ^ A((3) + -^(-3+10/3 2 +5/3 4 ) log + (1+0 Z ) 



1+/3 3 







4/3 : 



1-/3 2/3 2 



(3.16) 



with 



A{fi) = (1 + /3 2 ) 



4Li2 |iz^ +2Li2 f_lz^_31o * lo,^ 



1+/3 



1+/3 °l-/3 



2 log /3 log 



1+/3 
1-/3 



3/3 log 



1-/3 2 



4/3 log /3. (3.17) 



This may be compared with the corresponding formulas for the pair production 
of a colored fermion pair via vector and axial vector currents (the detail on the 
latter may be found in Appendix C) 



1 .... 33 + 22/32-7/34 1+/3 3(5 -3/3 2 ) 
= ^ A{J3) + ^3 (21 + 59/3 2 + 19/3 4 - 3/3 6 ) log i±| 



(3.18) 



(3.19) 



The dependence on /3 of the cross section normalized to its lowest order value 
is shown in Fig. 4. The corresponding quantities for fermion pair production 
are also shown for comparison. The correction for scalars is larger than that for 
fermions via vector current for all values of /3. The correction for fermions via 
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axial current is equal to that for vector current at the high energy limit (which 
reflects the chiral symmetry in this limit), but approaches the scalar curve at 
small (3. This latter behavior reflects the fact that spin does not play an important 
role in the nonrelativistic region. The difference of the vector current result is 
due to its S-wave dynamics contrast to the P wave behavior of the others. 

At the high energy limit (3 — > 1, one finds 

A = 3 , A V = A A = ^. (3.20) 

The correction for the scalar pair is four times larger. Near the threshold, one 
has 

2 2 

A~A A ~^-2, A v ~^-4. (3.21) 

The terms proportional to 1/(3 are the consequence of Coulomb gluon exchange. 

It is well known that perturbation expansion in a s breaks down very near the 
threshold and one has to sum over ladder Coulomb gluon diagrams [2]. For the 
case of scalar pair production with the P wave threshold behavior, the lowest- 
order cross section is propotional to /3 3 , which is modified to 1 by the 0(a s ) 
correction. The nonperturbative contribution further enhances the cross section 
and leads to a constant cross section at the threshold [7]. In the nonrelativistic 
approximation with the Coulomb potential (one gluon exchange approximation) , 
the cross section is 



which gives at the threshold 

n 2 C 2 a 2 s 

This effect becomes important only at very near threshold, (3 ~ 0.1, however. 



A = . (3.23) 
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If one approaches closer to the threshold, one will encounter P-wave (( bound 
states, similar to charmonium and bottomonium resonances. These states are 
expected to be less pronounced than the fermion bound states, since the resonance 
formation cross section is proportional to the derivative of the wave function at 
the origin and is suppressed by a factor of a 2 compared to the fermion bound 
states. In the Coulomb potential approximation, this subthreshold cross section 
is given by [7] 

<r± = ^f^HUC R * s ) 5 E ^S(s-M 2 ) , (3.24) 
LA n=2 n 

where M n = 2m — C^a 2 m/4n 2 . This Coulombic approximation ignores the effect 
of confinement and is not adequate especially for higher bound states. One needs 
to solve Schrodinger equation with a realistic color-force potential to obtain the 
cross section in the resonance region. 

If C, has a short lifetime with T <S a 2 m, it does not live long enough to form 
a bound state. There will be only a smooth bump or shoulder instead of the 
resonances at the threshold. The situation will be similar to the case for the top 
quark pair production [8] and the cross section can be calculated quite reliably 
within perturbative QCD. 



4. Jet configuration 

In this section we study the shape of the three-body final state ((g when 
the gluon is 'visible'. We compare various distributions with those for the more 
familiar final state of a fermion-pair plus a gluon. 

To define the final state configuration, it is convenient to use Lorentz-invariant 
kinematic variables. We use the following Dalitz variables 

2p-k _ 2p-k 
x = , x = , (4.1) 

s s 

for which the phase space density is constant. The ((g cross section is (we drop 
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the subscript ± from here on) 



da 

dxdx 



CrOs 1 
n (3 



2 



1 1\ l + (3 2 1 1-/3 2 
xx) 2 xx 4 



9 — 9 



(4.2) 

The infrared singularity arises when the gluon energy is very small, x ~ x ~ 0. 
At the high energy limit, additional singularity arises from the region of the 
phase space in which the gluon momentum is parallel to that of Q or Q. This 
collinear region is characterized as x ~ or x ~ 0. In these singular regions, it 
is practically impossible to observe the existence of the extra gluon in the final 
hadron system. We therefore treat the ^ + gluon system as a ( (f jet when the 
^ g invariant mass is sufficiently close to the ( mass. 

We thus divide the three-body phase space into the two-jet and three-jet 
regions. The two-jet region is defined as 



x < x r or x < x ( 



(4.3) 



and the rest is the three-jet region. The two-jet cross section is the sum of the 
CC cross section and the ((g cross section in the two-jet region. 

In the high energy limit j3 — > 1, the integrated three-jet cross section is found 
to be 



^3j/^o 



Crol s 



n 



2 Li 2 (x c ) + log 2 x c - ^- - 2(1 -x c ) log - — — + (3-2x c ) (1 -2x c ) 

D X r 



(4.4) 



whereas for fermions (both vector and axial) we have 



Crcx s 



n 



2 Li 2 (x c ) + log 2 x c - 



TV 



1(1 - x c )(l - 3x c ) log LJ?£ + hi - 2x c ) 
2 x c 4 



(4.5) 



In Fig. 5, we show the three-jet fraction as a function of x c - (Here and in the 
following figures, we plot normalized to the 0(a s ) total cross section, not 
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to the lowest-order cross section.) The three-jet fraction for the scalar is larger 
than that for the fermion at large x c but becomes smaller at small x c 0.1. The 
approximate form for small x c is 



z Crcx s 

O"3j/O"0 - 



n 



log x c + 2 log x c + 3 - — 

6 



(4.6) 



for scalars, and 

CRa s 



cr 3j /ao ~ 



TV 



,2 3, 5 ^ 

log X c + - log X c + - - — 



(4.7) 



for fermions. These predictions are not reliable at very small x c since higher 
order contributions become important. 

The three-jet fraction for finite m has a very complicated expression and we 
do not write down the analytic expression here. Instead, we show the numerical 
result for a^j a for several values of f3 in Fig. 6. For a fixed value of x c , the three- 
jet fraction rapidly decreases as (3 becomes smaller. Hard gluon emission is much 
suppressed when one approaches to the threshold. This suppression comes from 
the available phase space and from the structure of the gluon-emission vertex. 
The large QCD enhancement near the threshold can be attributed to two-jet 
configuration and one can assume that most events in the threshold region have 
no extra jets in experimental searches for new colored scalar particles. 

Toward the threshold region, the three-jet fraction in Fig. 6 for scalar pairs 
becomes similar to that for fermion pairs from axial current. The fermion pair 
from vector current shows a different behavior. This again reflects the fact that 
spin of the particle is not important in the nonrelativistic region. 

In jet analyses in e + e~ physics, one normally uses a jet-defining algorithm by 
combining momentum of hadrons and compares the result with the corresponding 
quantity in the quark-gluon system. A commonly adopted prescription uses 
the invariant mass of the hadron system to define a jet. To conform with this 
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algorithm, we may use the invariant mass variables 



(p + k) 2 _ (p + k) 2 

y = , y = 



(4.8) 



instead of x and x. These variables are related by 



m 2 m 2 

y = x-\ , y = x-\ 

s s 



The definition of three-jet events in terms of y is 



(4.9) 



y > y c and y > y c 



(4.10) 



The three-jet fraction as a function of y c is shown in Fig. 7. 

So far we have assumed that we can somehow distinguish a <^-jet and a gluon 
jet. If this is not the case, we need to treat the three final particles in a democratic 
way and define the three-jet region as 



y > y c and y > y c and y g > y c , 



(4.11) 



where y g = (p + p)/s. Since the region y g — > (or y g — > Am 2 / s) contains no 
singularity, this definition does not lead to significant modification of the result 
for small y c . In the high energy limit /?—>!, the three-jet fraction is found to be 



03j/crO 



7V 



2Li 2 



Vc 



i-yc 



+ log z 



i - y c ^ 



6 



Vc 



1 — 2v 

2(1 - 2y c ) log ^ + 3(1 - y c )(l - 3y c ) 

yc 



, (4.12) 



whereas for fermions (both vector and axial) we have 



03j/crO 



CROts 



TV 



2Li 2 



Vc 



i-yc 



+ log z 



i - y c _ 

y c 6 



|(1 - 2y c ) log + i(l - 3y c )(5 + 3y c ) 



. (4.13) 



15 



5. Conclusions 

We have studied 0(a s ) QCD correction to colored scalar pair production in 
e + e~ annihilation. Emphasis is put on gluon emission in the final state. We 
compared the configuration of the three-jet final state with those in fermion 
pair production and found that the two cross sections are quantitatively quite 
different. This has to be taken into account in the search and study of production 
of scalar particles such as squarks and leptoquarks at LEP-2 and future linear 
colliders. 

In calculating the helicity amplitudes e + e~ — > CCfi 1 ? we found that it is conve- 
nient to decompose the amplitude into two subprocesses e + e~ — > V* (V* = 7, Z) 
and V* — > ((g. The amplitudes take a simple form if one works in two different 
Lorentz frames, the e + e~ cm. frame for the former and the cm. frame in the 
latter. In particular, working in the "quasi- two-body" frame in the latter 1 — > 3 
subprocess allows one to take advantage of the strong constraints of rotational 
invariance on two-body states. 
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Appendix A. Helicity amplitudes for e + e — > ££g 

A number of formulations [9] to calculate helicity amplitudes [10] are found 
in the literature. Most of them are particularly suited for many-body final states 
and numerical evaluation of the amplitudes. The technique we employ here is a 
method [11] based on spherical- vector basis, which is suited for manual analytic 
calculations and gives amplitudes with clear physical interpretation. The method 
is most useful for 2 — > 2 body processes (also 1 — > 2, 2 — > 1). In these processes, 
angular momentum conservation strongly constrains the form of the amplitudes. 
In our technique, this constraint can be explicitly extracted in the form of a 
Wigner D function which depends on the scattering angles. This part represents 
the variation of the amplitude which is imposed by kinematics. The rest of the 
amplitude contains purely dynamical information. 

In general, the helicity amplitude for the process 

a(p a , A ) + b(pb, X b ) -> c(p c , A c ) + d(p d , X d ) 

(A n is the helicity of the particle a, etc.) in the cm. frame can be written as 

M = M(E c . m ., cos#; {A}) d J x °. Xf (e) e^^'V , (Al) 

where Aj = A a — Aft, A/ = A c — Ad, Jo = max(Aj, A/), and (0, <p) are the scattering 
angles (in the frame p a is along the z axis, p c is in the direction (0, </>)). The 
d xx , is the Wigner d function. The last two factors in this formula represent 
the "minimal" angular distribution imposed by the kinematics. The physical 
meaning of Jq is the smallest allowed angular momentum for the process. The 
cos6> dependence of M. comes from higher partial waves with J > Jo. If the 
process has only one partial wave J (e.g. for two-body decays or e + e~ — > 
the whole angular dependence can be extracted in the form of d J and the 0- 
dependent phase factor. 
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This technique can be applied to the process we consider if we note the 
following two points: 

1) The amplitude can be decomposed into two parts: (A) production of a 
virtual vector boson V*, e + e~ — > V* (V = 7 or Z); (B) the virtual vector 
boson decaying to ££g. The helicity amplitude for the whole process is the 
product of the two amplitudes with a fixed V* helicity (\y = ±1, 0), which 
are then summed over: 

M = -Y,—^—t E M(e-e+^V*)M(V*^((g). (A2) 
v s - m v A y =±i,o 

(Note that we need three polarization states for both 7* and Z*. A fourth 
(scalar) polarization state need not be included since V* couples to a con- 
served current.) 

2) The helicity amplitude for the second 1 — > 3 process V*(q) — > +((p) + 
g(k) can be treated in the same way as 2 — > 2 processes if one works in the 
"quasi- two-body" frame, the cm. frame of two of the final particles (e.g. 
CC) 5 not in the V* cm. frame. This frame is related to the e + e~ cm. frame 
by a boost along the gluon direction. 

First, we list the helicity amplitudes for V*(q, \y) — > ((p)+((p)+g(k, A) (note 
that the scalar particles have zero helicities). We take the V* (and g) direction 
as the z axis, and Q direction as (6>, 0). We follow the Jacob- Wick convention 
to fix the phase of the fermion/ vector wave functions. The amplitudes can be 
written as 

M(V*^((g) = g f g s T a M, (A3) 

with 

{Qce for V = 7, 

e(T 3? - Q( sin 2 6>w-) / sin 9 W cos 6\y for V = Z, 

and M. given by the following: 
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(A4) 



A n = Ay — A = ±2 



~ 2tv 2 sin 2 9 a\ j. , . 

(1 — r)(l — v 2 cos 2 6*) 



A n = Ay — A = ±1 



~ 2v / 2rt> 2 sin6>cos6> ,\ j, , . „. 

(1 — r)(l — v 2 cos 2 0) 



-^n = Ay — A = 



~ 2rw 2 sin 2 9 n . A . 

^ = "71 7n 2 2m "2. A7 

(1 — t)(1 — ir cos 2 6*) 

The helicity amplitude for e~e + — > V"* is simpler. We take the V* direction 
as the z axis and the e~ direction lying on the xz plane (with a positive a; 
component). The angle between the two directions is denoted by O. This choice 
is made to relate the two frames by the boost along the z axis, under which the 
helicity of V* does not change. The amplitude conserves electron chirality so 
that only A; = A(e _ ) — A(e + ) = ±1 is allowed. 

M = (-l) X Vg t V2Sdl iXv (G), 

with 

f-e forF = 7 , Ai = ±l, 

9i = i ( A8 ) 
I e(i> e =F fle) / sin 9\y cos 6>yj/ for V = Z, Aj = ±1. 

The explicit form of the <i functions needed is 

fi(l + A^cose) |A| = |a»| = 1, 
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Appendix B. Three-body phase space 



The differential cross section can be calculated from the helicity amplitudes 
with the Lorentz-invariant three-body phase space 



m = MVh -P-P-k) (2 $ 2p0 , (BD 



which in terms of the kinematic variables defined in Section 3 is 



SV (1 — T) 

d® 3 = 1 Q 247r 4 dr d cos 6 d cos d<f> , (B2) 



where 1 — /3 2 < r < 1, and u = (1 — 4m 2 /rs) 1 / 2 is the ( velocity in the CC c.m. 
frame. (We have integrated over one azimuthal angle on which the cross section 
has trivial dependence.) 

Two of the kinematic variables, r and 6>, specify the final state configuration 
and the other two determine the orientation of the configuration with respect to 
the initial beam axis. Integrating over G and 4> one has 



S?J (1 — T 1 

rf$ 3 = 1 > drdcosO . (B3) 
2567r J 



For analytic integration over the three-body phase space for m ^ 0, it is easier 
to use (cos 6>, v) as the integration variables, for which the integration region is 
< v < (3, -1 < cos0 < 1. 



^= s( lJ )v T-~4 ) dvdmse - <B4) 



The variables (x, x) have a special property that the phase space factor is 

20 



constant (Dalitz variables). 



cM>3 = ^dx dx . (B5) 

1287T 3 v ' 



The translation is made using 



r = l — x — x, vcos9 = —- — 3. (B6) 

X + X 



The physical phase space region is bounded by the inequality 



2 

xx(l - x - x) - — (x + x) 2 > . (B7) 
s 



Appendix C. Formulas for fermion pair production 



We collect results for the familiar process of fermion pair production e + e~ 
FF(g) [2,12,13] in this Appendix. These results have been independently calcu- 
lated by us. The FF7 and FFZ vertices in the lowest order have the form 



with 
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-ieQFlv and a ^—1h(vf - OF75) , (CI) 

cos uw sin 0\y 



v F = l - [T 3 (F L ) + T 3 (F R )} - Q F sm 2 9 w , (C2) 
a F = l[T 3 (F L )-T 3 (F R )] . (C3) 



At 0(ct s ), the correction factor is common for the photon and the vector part of 
the Z current. The axial part of the Z current is modified differently if rap 7^ 0. 
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We write the general Lorentz structure of the vector and axial vertices for on-shell 
fermions 



r M = i^V)7M + F 2 v (q 2 )ia^/m F 
for the vector current and 



= FiVh/iTS + F 2 A (q z )q^ 5 /m F 



At 2\ 



(C4) 



(C5) 



for the axial current. We have retained only terms appearing at 0(a s ), i.e., 
terms allowed by CP invariance and vector current conservation. We normalize 
the vertices such that F^ ,A = 1 at the lowest order. F 2 V,A do not appear in this 
order. 



We write the 0{a s ) corrected vertices as 



(C6) 
(C7) 



We find the renormalized vector form factors for q 2 > 4m F 



1 + /3 2 . >n 
log - — - + 1 ) lo 



+ 



+ 



2(3 

1+/3 2 

(3 

1 + 2/3 2 



Lio 



1-/3 ' 7 "° A 2 

1+/3 l+(3 1 2 l+(3 tt 2 
+ log log log 

1+/3 J S 2(3 S 1-/3 4 S 1-/3 3 



log - — - - 2 



1-/3 



+ Z7T 



1 + /3 2 



2/3 



log 



777-1 



A 2 



+ log 



4/3 2 



1 + 2/3 2 
2/3 



J2 — 



1-/3 2 



log + Z7T 



(C8) 



(C9) 



4/3 V. 1-/3 

where /3 = (1 — Am 2 F / q 2 ) 1 ^ 2 . We have used on-mass-shell renormalization con- 
dition such that F^(0) = 1. Once the prescription for the vector vertex is fixed, 
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there is no freedom to choose a renormalization condition for the axial vertex be- 
cause of electroweak gauge symmetry. In particular, F^(0) deviates from unity. 
The axial form factors are found to be 



ft 



log - — - + 1 I log 



+ 



+ 



2/3 
1 + /3 2 



1-/3 



A 2 



/3 

2 + /3 2 l+p 



, . . 1-/3 \ , 1+/3 1 1+/3 1. 2 1+/3 7T 2 

Lio — - + log — - log — log — - H 

" 1 1+/3; 8 2/3 S 1-/3 4 8 1-/3 3 



2/3 



log 



1-/3 



-2 



+ Z7T 



^6og^ + logi^V^ 

2/3 V A 2 2/3 



M _ 

J2 — 



1-/3 
4/3 



2 r 



(2 + /3^)(-logl±| + Z7r)+2/3 



(CIO) 
(Cll) 



The second form factor represents the pseudoscalar part of the axial current 
and does not contribute to the reaction we are interested in as long as the electron 
mass is neglected. 

The differential cross section for the 3-body FFg final state is for the vector 
current 



dal 



dr d cos 9 



v C R a s v(l - t) 



x 



/3(3-/3 2 ) 
l + v 2 cos 2 9 ^ 2(3-/3 2 )™ 2 sin 2 # 



l-^cos 2 ^ (l-r) 2 (l-?; 2 cos 2 6>) 2 



(C12) 



and for the axial current 



dot 



dr d cos 9 



A C R a s v(l - r) 
2 - /3 2 + v 2 cos 2 9 



x 



4/3 2 ri> 2 sin z 6> 



+ 



l-v 2 cos 2 9 (1-t) 2 (1-v 2 cos 2 9) 2 



(C13) 
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with the lowest order cross sections 



v 8ird R a 2 u2 (3(3 -(3 2 ) 



3s 



l v±- 



3s 



Here Hy± and Ha± are given by 



Ha± 



-Qf + 



(v e Ta e )v F s 



cos 2 9w sin 2 9w s ~ TO | 

(v e Ta e )aF s 



cos 2 9w sin $w s - m z 



(C14) 
(C15) 



(C16) 
(C17) 



In terms of the Dalitz variables, we have for the vector current 



da\ 
dx dx 



V C R a s 1 
a0± — p 



3-(3 2 



1 1 
- + - 

x x 



1 + (3 2 1 1 - f3 2 ( 1 1 



2 xx 



and for the axial vector current 



do± _ _ A C R a s 1_ 
n 13 



dx dx 



3-(3 2 (x x\ l-[3 2 
~j0 2 ~\x Jr x~ 1 1 



1 1 

- + - I + 

X X 



2(3 2 

1 + /3 2 1 l-[3 2 
2 xx 4 



2 



1 1 

/j^2 rjp 2 



(C18) 



(C19) 



The total 0(a s ) correction can be written as the sum of the three contribu- 
tions, virtual, soft, and hard corrections. We list each correction term for the 
scalar-pair production for comparison 



A v i rtu al = -pA v {p) + log — jj - 2 



(C20) 
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^soft 



Ahard = 



^(/3)-^3(3 + /3 2 )(l-/3 2 )log 1 + /3 



(C21) 



+ -^(3 + 7/3 2 ), (C22) 



1-/3 2/3 2 



where we have collected the "dilogarithmic" part (dilogarithm and double log 
terms) into the functions Ai (i = v, s, h) 



A v (f3) = 



MP) 



A h (P) 



1+/3 2 , 1+0 



2 l0g T^ + /3 J l0g A^ 



m 



+ (l + /3 2 ) 
/1+/3 2 1+P 



\l+p) 1+P 1-/3 4 s 1-/3 3 



(C23) 



V. 2 ° 1-/3 
+ (l + /3 2 ) 



lo S i — 3 _ P lo § 



4^ 
A 2 



/l-/3\ 2/3 1+/3 1 2 l+/3 tt 2 

Ll2 il^J- l0g l^ l0g I^-4 l0g 1^/3 ~ T 



(C24) 



1+/3 2 , 1+/3 



log 



1-/3 



/3 log 



m 
Aj 2 



+ (l+/3 2 ) 



2 02 (irf ) + 2 02 



1, 2 1+P ^ 

+ o lQ g 



1-/3 6 



1-/3 
1+/3 



3/3 log 



1-/3 2 



i 2 i !+/3 
log log 

S 1+/3 S 1-/3 



4/3 log /3. (C25) 



The function A{(3) in (3.17) is the sum of these three 



A{(3) = A v {(3) + MP) + MP) 



(C26) 



Turning to the fermion-pair production, it is found that the dilogarithmic 
terms are exactly the same as for the scalar-pair production. For the vector part 
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we find 



^virtual = **fY + J^2^f¥ 



-A (3) I ffl^ wl + g 



Ar oft =^) + iio g i±|, 



/3 ftV ^ 8/3(3 - /3 2 ) 4(3-/?" 
and for the axial part 



A4 tu , 1= ^ (/3 ) + ^.o g i±|-2, 
A^ = ^. (/3 ) + ilogi±|, 



A^ ard = ^fcGS) + 3^3 (21 - 5/3 2 + 3/? 4 - 3/3 b ) log - 
+ I ^2(-21 + 62/3 2 + 3/3 4 ). 
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Figure Captions 

Fig. 1. Feynman diagram for e + e~ — > (( a t the lowest order. 
Fig. 2. Feynman diagrams for e + e~ — > CC at C(q; s ). 
Fig. 3. Feynman diagrams for e + e~ — > (Xg. 

Fig. 4. Total C(a s ) correction to the pair production cross section of scalar pair 
(solid); fermion pair via vector current (dash); fermion pair via axial 
current (dotted). All curves are for particles in the fundamental color 
representation (Cr = 4/3) and the strong coupling constant a s = 0.12. 

Fig. 5. Three-jet fraction for scalar (solid) and fermion (dashed) at the high 
energy limit, for a s = 0.12 and Cr = 4/3. 

Fig. 6. Three-jet fraction for scalar (solid), fermion- vector (dashed), and fermion- 
axial (dotted) for (3 = 0.4, 0.6, 0.8, and 1.0, a s = 0.12 and Cr = 4/3. 

Fig. 7. Same as Fig. 6, but as a function of a jet-defining variable y c . 
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